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Abstract 

In this paper it is shown that Tolman's law can be derived from rel- 
ativistic kinetic theory applied to a simple fluid in a BGK-like approx- 
CO. imation. Using this framework, it becomes clear that the contribution 

of the gravitational field can be viewed as a cross effect that resembles 
the so-called Thomson effect in irreversible thermodynamics. A proper 
generalization of Tolman's law in an inhomogeneous medium is formally 
C ' established based on these grounds. 
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^ ■ 1 Introduction 

In 1930 Richard C. Tolman showed that thermal equilibrium can persist within 
a temperature gradient provided that a gravitational field is present, satisfying 
the relation: 

L*L =9 (i) 

T dr c 2 

which is often referred to as Tolman's law Q]. Later, in 1940, Eckart proposed 
that the heat flux in a special relativistic fluid should be proportional to the 
temperature gradient and to the hydrodynamic acceleration of the fluid [5] ■ The 
corresponding constitutive equation, derived from phenomenological arguments 
reads 

T 



Jfa = -h^Kr { T v + -a, ) (2) 
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where h^ v is the spatial projector given by 

bT = gT + \WU V (3) 
& 

with W L being the hydrodynamic four- velocity, U^U^ = — c 2 and kt is the ther- 
mal conductivity. It seems attractive to relate Eckart's constitutive equation 
to Tohnan's law since requiring a vanishing heat flux for equilibrium leads to 

ldT = _ae 

Tdx 1 c 2 k ' 

which strongly resembles Eq. ([1]). However, in the presence of stresses the 
hydrodynamic (macroscopic) acceleration is not equal to gravity. Moreover, 
the coupling expressed by Eq. @ does not relate a thermodynamic flux with a 
canonical force, since a v is not the spatial gradient of a state variable as required 
by irreversible thermodynamics. Also, this coupling has been proven to be the 
cause of generic instabilities, first identified by Hiscock and Lindblom in 1985 

mm- 

On the other hand, it is well known that kinetic theory provides a framework 
in which constitutive equations can be established from first principles, in com- 
plete consistency with linear irreversible thermodynamics in the Navier-Stokes 
regime. In this paper, it is shown that a non-vanishing contribution to the heat 
flux due to a classical gravitational field arises within the framework of relativis- 
tic kinetic theory. Such effect strongly resembles the Thomson effect present in 
plasmas under the influence of an electrostatic field. It becomes equivalent to 
Tolman's law for a homogeneous medium at mild temperatures in the presence 
of a gravitational field. 

In section 2 the basic principles of general relativistic kinetic theory within 
the BGK approximation are described. Section 3 is devoted to the solution 
of Boltzmann's equation and to the formal expression for the heat flux, in the 
presence of a gravitational field, using the Chapman-Enskog expansion to first 
order in the gradients. In section 4, the corresponding transport coefficient 
is calculated showing that Tolman's law follows by requiring a vanishing heat 
flux in a homogeneous medium. A proper generalization of Tolman's law in 
an inhomogeneous medium is also formally established. Conclusions and final 
remarks are included in section 5. 



2 Relativistic kinetic theory 

The standard method for establishing transport equations as well as constitutive 
relations for dissipative fluxes is well known, and has been thoroughly discussed 
by several authors. The formalism has been successful in many applications, 
including relativistic fluids [5J [5J H] • In such program a collisional term for the 
kinetic equation must be specified as a starting point. In a relaxation time 
approximation, using the so-called BGK collision kernel for the relativistic gas 
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proposed by Marie [5], the kinetic equation reads 

v a la + a»^- = -t— (5) 

where / is the single particle distribution function, /(°) the local equilibrium 
(Juttner) distribution [5] [5], i> M is the molecular four- velocity and r is a pa- 
rameter of the order of the collisional time. Since molecules follow geodetic 
trajectories, the molecular acceleration can be written in terms of the cur- 
vature induced by a gravitational field as 

= -r> a «" (6) 

Thus, Eq. §5§ becomes 

.at _ r M ..a v d f _ / ~ / ° ( 7 \ 
V T,a L af }V V ^ - ^ , (I) 

According to the Chapman-Enskog method, it is possible to establish a formal 
solution to the Boltzmann equation in the Navier-Stokes regime that reads: 

/ = /«>)+ /(i) (g) 

where is a correction to local equilibrium to first order in the gradients. 
Substituting expression (|5J) in Eq. ([?]). and retaining only first order terms in 
the gradients leads to 



/P) = -r{«-A»'-V^| (9) 



Q/3 " dv^ J 

Using the local equilibrium hypothesis, we can write the derivatives of f(°> in 
Eq. © as: 

/(0) = ^ T + d J^n + (10) 

where T is the local temperature and n is the local number density. Recalling 
that the Juttner function reads: 

f (o) = « exp ( fU ) 

a direct calculation leads to 



r(0) 



f(o) J IV 




T 



(12) 



where A'„ 



is the nth modified Bessel function of the second kind, Q 



K3 ^ f j and 7 = is the standard Lorentz factor for the molecules' chaotic 



K, 



velocity k (for more details see Ref.[7]) 
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Since the molecular four-velocity that appears in Eqs. (|lltil2p can be written 

as 



the first term of Eq. © can be expressed as 



(13) 
(14) 



or 



v a f. 



(0) 




where §i\ CF indicates a derivative calculated in the comoving frame, where 
locally U v = [0, c]. Following Hilbert's standard procedure [5], Euler equations 
are now used in order to express the time derivatives in Eq. (|15p in terms of 
the gradients. We recall that Euler's equations in the comoving frame read (see 
Ref. 0): 



dU l 
dt 



CF 



dn 






~dt 


CF 


T 




n, v 


— H 

T 
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dT 




~dt 


CF 



-c 2 r 
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dU A 



dt 



= 



CF 







(16) 
(17) 
(18) 



In Eq. (fTT|) use has been made of the fact that all four components of the 
projector vanish as well as r| 4 in a static metric. 

The last term of Eq. ([9]) is easily established, and for a newtonian metric 
reads: 



t a/(°) 



We are now in position to address the main task of this work namely, to analyze 
the heat flux associated to a gravitational field applied to a simple fluid. 



(19) 
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3 Heat flux 



In this section attention will only be paid to the gravitational terms present 
in the first order in the gradients correction to the equilibrium distribution 
function. Introducing Eqs. (P~5|) to (fH))) in Eq. © and ignoring all terms 
related to non-gravitational potential gradients one obtains 

f$ = -^f i0) { (-r&) + t,^)} , (20) 

where p = + ^ . This expression clearly vanishes in the non-relativistic limit 
and thus no heat flux is induced due to a gravitational field in a simple non- 
relativistic fluid, as expected. However, in the relativistic case and using the 
newtonian approximation for the metric, the following expression can be written: 

and thus 

fCO rti°)~ W J 7(fc) A / on 

4] = - TfU ^l^\g^)- 1 j (21) 

We now recall that the local heat flux in a relativistic fluid is expressed in 
the Navier-Stokes regime as [5]: 

j( Q] =mc*J k e ^h 2 (k) d*K (22) 

so that the thermal dissipation term that arises from a gradient in the gravita- 
tional potential can be identified as: 

43..] = -"V 1 / f (0) ^h\ {^fj - l} W (23) 

Notice that the determinant of the metric is not included as a factor in the 
integral since within this approximation deviations from a unitary value are 
second order and higher in the gravitational potential gradient. The invariant 
velocity volume element d*K is established in Refs.|10[ ITT] . Details on the 
procedure to evaluate this type of integrals can also be found in these references. 
The constitutive equation for the gravitational term in the heat flux can then 
be written as 

J [Q, 9 ] = - W (24) 

where 

£ -=^{ fc+ srir s (3)}' (25) 



Figure 1 : Dependence with z of the coefficient that couples the heat flux with the 
gravitational field. The coefficient tends to zero in the small z, non-relativistic, 
limit and asymptotes to a constant value for z>l (ultra-relativistic systems). 




The behavior of Lt q can be observed in Fig. 1 where f g (z) = Lx g / (kTrn) is 
plotted as a function of z. 

To the authors' knowledge this calculation is novel. In the next section we 
shall analyze the low temperature limit of Eq. (|24|) and relate the corresponding 
result to Tolman's law. 



4 Transport coefficients and Tolman's law 

For z€l, the expansion in Taylor series of the expression in brackets in Eq. 
([2"5|l reads: 

i / 1 \ n 

(26) 



5z 
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and to lowest order in z one obtains 



\\ )-§(.'-=■ + ...) 



T t 5 nk 2 T 2 , 

If the fluid is assumed to be homogeneous (constant density), the heat flux is 
expressed as follows: 



IQ] 



rp,e 

-L T t— - L Tg g 



(28) 



where it is easily seen that since Lt 9 scales as z, the cross effect that couples 
heat with gravity is only present in the relativistic case. 
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In order to relate this result with Tolman's law we assume a vanishing heat 
flux as a requirement for equilibrium so that 

^ = -^V (29) 

1 LtTT 

Since the heat conductivity coefficient to lowest order in z is given by [5] 

5 nk 2 T 2 

Ltt = t; t, (30) 

2 m 

one finally obtains the standard form of Tolman's law: 

T' 1 1 

- = ^ (31) 
It is worth to notice that for a non-homogeneous medium, the heat flux 
expression involves three thermodynamic forces 

e T' 1 w l 

J[Q] = -L TT — — LnT ~ ~ L Tg9 ( 32 ) 

Therefore in the general case of a simple relativistic non-degenerate and inho- 
mogeneous fluid, the following generalization of Tolman's law is obtained: 

T' e n- 1 

Lrgg = —Ltt-^ — L n r — (33) 
in 

The complete expressions for the transport coefficients Ltt and L„t have al- 
ready been calculated in reference [5] using the BGK approximation for the 
kernel in Boltzmann's equation. Namely, 

nk 2 T 2 r 




(34) 



-^{^-(-f§)>-©)}. » 

It is important to point out that, in the case of an inhomogeneous, non- 
isothermal fluid in the presence of a gravitational field, the contribution for the 
total heat flux is of order zero in z for the temperature gradient and of first 
order in z for the density and gravitational potential gradients. Therefore, the 
temperature gradient makes the most important contribution in the heat flux 
generation. 
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5 Discussion and final remarks 



It is well known that a heat flux arises when a temperature gradient is present 
in a simple fluid. Richard C. Tolman showed, based only on phenomenological 
arguments, that a gravitational potential can balance a temperature gradient 
in order to achieve thermal equilibrium. In this paper it has been shown, using 
relativistic kinetic theory arguments that, in a relativistic inhomogenous simple 
fluid, a heat flux can be generated not only by a temperature and/or density 
gradient but also by a gravitational field. This is expressed by Eq. ([32| and 
the corresponding gravitational transport coefficient has been established in the 
Navier-Stokes regime using the first order Chapman-Enskog method to solve 
Boltzmann's equation. 

A further analysis of expression (|3"2"j) . illustrates the fact that the heat flux 
is only related to canonical forces. It is important to notice that, in a continu- 
ous medium, the hydrodynamic acceleration is not equal to the single particle 
acceleration. In general Tolman's law should not be identified with Eq. (j4|. 

The expression for the total heat flux shows that two cross effects are present 
in a single component relativistic fluid. This feature is completely absent in the 
single component non-relativistic treatment in the presence of external fields. 
The present calculation suggests that a non-vanishing diffusive flux may also 
exist in the relativistic regime implying symmetry (Onsager's) relations for the 
transport coefficients. 

Kinetic theory provides a very powerful tool to explore the transport prop- 
erties of gases in the presence of strong gravitational field barely studied in the 
past. The extension of this result to non-static space-times is a very interesting 
subject which will be addressed in the near future. 
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